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Summary. — I discuss a few selected topics related to the reconstruction of the 
mass of the top quark at hadron colliders. In particular, the relation between the 
measured top mass and theoretical definitions, such as the pole or MS mass, is 
debated. I will also summarize recent studies on the Monte Carlo uncertainty due 
to the fragmentation of bottom quarks in top decays. 

PACS 14 . 65 . Ha - Top quarks. 

PACS 12.38.Bx - Perturbative calculations. 



The top quark is an essential ingredient of the Standard Model of the fundamental 
interactions and the determination of its properties provides important tests of the strong 
and clcctrowcak interactions. In particular, the top quark mass (m*) plays a crucial role 
in the precision tests of the Standard Model, as the combination of m t with the W-boson 
mass mw constrains the mass of the Standard Model Higgs boson (see, e.g., the updated 
results from the LEP electroweak working group [1]). Moreover, the top mass, along 
with mw and the mass of a possibly discovered Higgs boson, is a relevant parameter 
to discriminate among New Physics scenarios, such as the Minimal Supcrsymmetric 
Standard Model [2]. The latest top mass measurement at the Tevatron accelerator is 
m t = (173.3 ± 1.1) GcV [3]. However, it is often unclear how to relate the measured m t 
to a consistent mass definition. 

In fact, from the theoretical point of view, several mass definitions exist, according to 
how one subtracts the ultraviolet divergences in the renormalized self-energy S(p, m t , n), 
where p is the momentum of the top quark, assumed to be off-shell, and fi the renormal- 
ization scale. The pole mass, often used in decay processes of on-shell particles, is thus 
defined by means of the conditions E(p) = and dl^/dyf — for rf= mt- The pole mass 
works well for leptons, such as electrons, whereas, when dealing with heavy quarks, it 
presents the so-called non-perturbative ambiguity, i.e. an uncertainty Am ~ Aqcd- In 
other words, after including higher-order corrections, the self-energy, expressed in terms 
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of the pole mass mt, exhibits in the infrared regime a renormalon-like behaviour: 

(1) E(mt) ~ m* a s +1 (W n! ' 

n 

where /3q is the first coefficient of the QCD /3-function. 

Another definition of top mass which is often used is the MS mass, namely m t (/K), 
corresponding, e.g. in dimensional regularization with D = 4 — 2e dimensions, to sub- 
tracting off the rcnormalizcd self-energy the quantity 1/e + "/ F — ln(47r), "f F being the 
Euler constant. The MS mass is a suitable mass definition in processes where the top 
quarks are off-shell; in fact, by means of this definition, one is able to reabsorb in mt(/z) 
contributions ~ ln(/x 2 / m 2 ) , which are large if \i is taken of the order of the hard scale 
Q and Q ^> mt- For ti production at threshold, nonetheless, the MS mass is not an 
adequate mass scheme, since it exhibits contributions ~ (ag/v) , v being the top quark 
velocity, which are enhanced for v — > 0. Of course, one can always relate pole and MS 
masses by means of a relation like [5] : 

(2) m t = m t (/i) [l + as(n)ci + a 2 s {ji)c 2 + ...], 

where the coefficients Cj depend on ln[m 2 (/x)//x 2 ]. 

Besides pole and MS schemes, other mass definitions for heavy quarks have been 
proposed. For example, in the potential-subtracted (PS) mass a counterterm Smt is 
constructed in such a way to remove the renormalon ambiguity in the infrared regime at 
the factorization scale \xf- It is given by [4]: 

(3) rn ttPS (iJ. F ) = m t - Sm t (p F ), 
where the subtracted term reads 

1 f d 3 q ~ 



(4) 5m( PF ) - - / 7^3^(9), 

2 J\q\<flF ( 27r ) 

V(q) being the Fourier transform of the ti Coulomb potential. More generally, one can 
define several short-distance masses in terms of a parameter, usually called R, corre- 
sponding, e.g., to the factorization scale /ip or the MS mass, so that the pole mass mt 
can be expressed in terms of a independent renormalized mass and a counterterm [6] : 

(5) mt = mt(R,(J.) + Smt(R,fJ,). 

In Eq. (5) Smt(R, n) can be expanded as a series of the strong coupling constant, whereas 
m t {R, fi) satisfies a renormalization-group equation: 

(6) 5m t (R,n) = R^Y^ankOsifi) , = -^hMl 

with 7[ag(/x)] playing the role of the anomalous dimension [6]. 

In order to make a statement on the top quark mass which is extracted from the 
data, in principle one should compare the measurement of a quantity depending on the 
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Fig. 1.: Left: Extraction of the top pole mass by comparing the DO cross section with 
the computations [8, 9]. Right: Determination of the MS mass following [9]. 



top mass with a calculation which uses a consistent definition, e.g. the pole or the MS 
mass. The DO Collaboration [7] compared the measured total ti cross section with the 
calculations [8] and [9]. Ref. [8] computes the ti cross section at next-to-leading order 
(NLO) and resums soft /collinear contributions in the next-to- next-to-leading logarithmic 
(NLL) approximation. The calculation [9] is NLO, resums the Sudakov logarithms of the 
top velocity a| ln fc v, with k < 4, and includes Coulomb contributions ~ 1/v and ~ 1/v 2 
to next-to- next-to- leading order (NNLO). Comparing the experimental cross section with 
the formulas in Rcfs. [8, 9], one can extract the top mass employed in such calculations. 
Using the pole mass, which is a reasonable assumption since top quarks are slightly above 
threshold at the Tevatron, one obtains the following results [7]: m t = 171.5^g'g GeV, 
according to [8], m t — 173.1^g'g GeV when using [9]. 

Furthermore, in Ref. [9] the computation was carried out in the MS renormalization 
scheme, and fh t (rht) was extracted from the comparison with the measured cross section. 
Fig. 1, taken from Refs. [7, 9] shows the main results of such a comparison. From fh t (rh t ) 
one can determine the pole mass mt : the central values are reported in Table I, whereas 
the errors are quoted in [9] . Ref. [9] also found out that using the MS top mass leads to a 
milder dependence of the cross section on factorization and renormalization scales, with 
respect to the pole mass. This result is indeed quite cumbersome, since, in principle, 
at the Tevatron top quarks are almost at threshold; the possible full inclusion of NNLO 
terms should probably shed light on the uncertainty when using pole or MS masses. 





m t {m t ) 






m t 


1 L0 1 


159.2 GeV | 


159.2 GeV | 


| NLO 


159.8 GeV | 


165.8 GeV | 


| NLO+ NNLO approx. | 


160.0 GeV | 


168.2 GeV | 



Table I.: Extracted values of the pole and MS top mass comparing the computation [9] 
with the DO ti cross section. See [9] for the errors on the quoted numbers. 
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Although the DO determination of the top mass by using fixed-order and possibly 
rcsummed calculations is surely very interesting, the standard Tcvatron analyses [3], 
such as the template or matrix-element methods, are driven by Monte Carlo parton 
shower generators, such as HERWIG [10] and PYTHIA [11]. Such algorithms simulate 
multiple radiation in the soft/collinear approximation and are possibly supplemented 
by matrix-element corrections to include hard and large-angle emissions. Moreover, the 
yielded total cross section is LO, whereas differential distributions are equivalent to 
leading- logarithmic soft/collinear approximation, with the inclusion of some NLLs [12]. 
The hadronization transition is finally implemented according to the cluster model [13] 
in HERWIG and the string model [14] in PYTHIA. 

The template and matrix-element methods rely on the Monte Carlo description of top 
decays. As discussed above, when looking at top decays near threshold, the reconstructed 
top mass should be close to the pole mass and the world average actually agrees, within 
the error ranges, with the pole mass extracted from NLO cross section calculations. 
However, as parton shower generators are not NLO computations, there are several 
uncertainties which affect the correspondence between the mass which is implemented 
and any mass definition and renormalization scheme. Monte Carlo algorithms neglect 
width and interference effects, but just factorize top production and decay; this is a 
reasonable approximation as long as one sets cuts on the transverse energy of final- 
state jets much higher than the top width T t . Nevertheless, the quoted systematic and 
statistical errors on the top mass [3] are competitive with the top width. Furthermore, 
top decays (t — > bW) in both HERWIG and PYTHIA are matched to the exact tree-level 
t — > bWg calculation [15, 16], but the virtual corrections in the top-quark self-energy, 
which one must calculate to consistently define a renormalization scheme, are included 
only in the soft/collinear approximation by means of the Sudakov form factor. Also, 
other uncertainties are associated with the flow of the colour of the top quarks and with 
the hadronization corrections, since all measured observables are at hadron-level. 

In principle, higher-order calculations are available even for top decays: Refs. [17, 18] 
calculated the 6-quark spectrum in top decay at NLO and included NLL soft/collinear 
resummation in the framework of perturbative fragmentation functions. Nevertheless, 
such calculations are too inclusive to be directly used in the Tevatron analyses; also, 
the results are expressed in terms of the 6-quark (£>-hadron) energy fraction in top rest 
frame, which is a difficult observable to measure. Ref. [19] recently computed several 
quantities relying on top decays at NLO using the top pole mass: a comparison of this 
calculation with Tevatron or LHC data, though not carried out so far, should possibly 
provide another consistent determination of the pole mass. 

An attempt to associate the Monte Carlo top mass with a formal mass definition has 
been carried out in [6, 20, 21, 22] in the framework of Soft Collinear Effective Theories 
(SCET). According to this approach, valid in the regime Q > m ( > T f > Aqcd, where 
Q is the process hard scale, one can factorize the double-differential cross section in terms 
of the top (Mt) and antitop (Aff) invariant-mass squared as follows: 

(7) dM fdM? =a ° H Q Hm J d£ + drB + (e+,T t )B^(£-,r t )s(i + ,e-). 

In Eq. (7), Hq and H m are hard-scattering coefficient functions at scales Q and mt, 
B± arc the so-called heavy-quark jet functions, describing the evolution of top quarks 
into jets, S(£ + ,£~ , /j,) is named soft function and is a non-perturbative fragmentation 
function, depending on soft radiation and ruling the dijet and mass distributions. When 
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writing a factorization formula like in (7), large logarithms of the scale ratios, such as 
ln(Q/m t ) or ln(mt/r t ), clearly arise: the jet function has been resummed to NNLL 
in e + e~ annihilation, the soft one to NLL [20]. The peak value of the distribution 
d 2 a/(dM^dM?), displayed in Fig. 2 in the NLL approximation for e + e~ — > tt processes, 
is independent of the mass scheme and can be expressed in terms of a short-distance 
mass to, r t , Q and Aqcd: 

(8) Mf cak = m + r t ( Cl « s + c 2 al + ...) + 

TO 

where the term ~ T t depends on the jet function and the one ~ Aqcd on the soft 
function. The jet mass is thus constructed as a short-distance mass, defined following 
Eq. (5), with the parameter R ~ F t and the counterterm 8m j depending on the jet 
function. Expressing M t pcak in terms of the pole (m t ) and jet (toj) masses, one can 
finally relate mt and toj(^) [21]: 

(9) m t = mj M + e^T t ^^ (in £ + |) + 0{a%). 

One can note that a correction ~ T t arises in the difference between the pole and jet 
masses: in fact, it was pointed out above that the top width was one of the ambiguities 
when associating the top mass reconstructed from final-state decay obscrvables, e.g. the 
jet mass, with the pole mass. Rcf. [22] then assumed that the jet mass should correspond 
to the mass in the event generators and measured at the Tevatron with [i being of the 
order of the hadronization scale (shower cutoff), i.e. [i ~ Qo ~ 1 GeV. Using the R- 
cvolution equations (5), with toj(Qo) — 172 GeV, one can finally obtain the MS mass 
m t (m t ) ~ 163.0 GeV [22] (see Fig. 2). Refs. [6, 21, 22] seem therefore to indicate a 
possible strategy to relate the mass parameter in the Monte Carlo codes to pole and 
MS masses. However, as discussed also in [22], the situation at hadron colliders is more 
complicated and, before drawing a final conclusion, one would need to take into account 
effects, such as initial-state radiation, hadronization and underlying event, which are 
included in Monte Carlo event generators, but not yet in the factorization formula (7), 
which has been worked out for e + e~ annihilation. Further studies aiming at extending 
the factorization (7) to hadron colliders are certainly worth to be pursued. 




Fig. 2.: Left: Double differential resummed cross section with respect to the top and 
antitop invariant masses. Right: i?-evolution from the Tevatron top mass, corresponding 
to R ~ 1 GeV, to the MS mass, i.e. R ~ m t (rht). 
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In fact, the theoretical error on the top mass determination is a crucial issue, even 
for the sake of the Standard Model precision tests. Ref. [23] presents an analysis aimed 
at assessing the uncertainty due to the treatment of bottom-quark fragmentation in top 
decays in HERWIG and PYTHIA. For this scope, one should tune the cluster and string 
models to the same data sets: Ref. [24] found out that the default parametrizations are 
unable to reproduce LEP and SLD data on i?-hadron production in e + e~ annihilation, 
but it was necessary to fit the hadronization models to obtain an acceptable description of 
such data. As in [23], an interesting observable is the invariant mass mse in the dilepton 
channel, where B is a ^-flavoured hadron in top decay and £ a lepton from W decay 
W — > tu. A reliable description of mst is essential in the study [25] (see also Chapter 8 
in Ref. [26]), wherein the top mass is reconstructed by means of the mj/^i spectrum, 
with the J/tp coming from B — > J/ip decays, and in the analysis [27], which investigates 
semileptonic decays B — > fiX and extracts the top mass from the m^i distribution. 
Moreover, the i?-lepton invariant mass was employed in [28] to gauge the impact of the 
implementation of hard and large-angle radiation in the simulation of top decays; also, 
rriBt is one of the observables calculated in the NLO approximation in Ref. [19]. 

Fig. 3 presents the vast spectra using tuned and default versions of HERWIG and 
PYTHIA, whereas Tables II and III quote the first few Mellin moments of vtibi for 171 
GcV < m t < 179 GeV, according to tuned HERWIG and PYTHIA. From Fig. 3 one 



| m t (GeV) || 


{mat) (GeV) | 


(m 2 m ) (GeV 2 ) | 


{ml,) (GeV 3 ) | 


(m 4 m ) (GeV 4 ) | 


| 171 


78.39 


7.01 x 10 3 


6.82 x 10 s 


7.02 x 10 8 


| 173 


79.52 


7.22 x 10 3 


7.12 x 10 s 


7.43 x 10 s 


| 175 


80.82 


7.45 x 10 3 


7.46 x 10 s 


7.91 x 10 8 


| 177 


82.02 


7.67 x 10 3 


7.79 x 10 5 


8.37 x 10 8 


| 179 


83.21 


7.89 x 10 3 


8.13 x 10 5 


8.86 x 10 8 



Table II.: First four Mellin moments of the rriBi spectrum according to HERWIG tuned 
to LEP and SLD data, as in [24], for 171 GcV< m t < 179 GeV. 
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| m t (GeV) || 


(m Be ) (GeV) | 


(m% t ) (GeV 2 ) | 


(m% e ) (GeV 3 ) | 


(m*Bt) (GeV 4 ) | 


| 171 


77.17 


6.85 x 10 3 


6.62 x 10 5 


6.81 x 10 s 


| 173 


78.37 


7.06 x 10 3 


6.94 x 10 5 


7.23 x 10 8 


| 175 


79.55 


7.27 x 10 3 


7.25 X 10 s 


7.67 x 10 8 


| 177 


80.70 


7.48 x 10 3 


7.56 x 10 5 


8.12 x 10 s 


| 179 


81.93 


7.71 x 10 3 


7.91 x 10 5 


8.61 x 10 8 



Table III.: As in Table II, but employing the PYTHIA code. 



Solid: HEHWIG 



Dots: PYTHIA 



Fig. 4.: Average value of the B£ invariant-mass spectrum in terms of nit, according to 
the best-fit straight lines yielded by the tuned versions of HERWIG and PYTHIA. 



learns that the tuning has a relevant impact on such spectra; Tables II and III show that 
visible differences between HERWIG and PYTHIA are still present, even after the fits to 
LEP and SLD data. To relate such a discrepancy to an uncertainty on the Monte Carlo 
top mass, one can try to express (msi), or even the higher moments, in terms of m t by 
means of a linear fit using the least-square methods. The best-fit straight lines read: 

(10) (™>bi)-r — —25.31 GeV + 0.61 m t ; 8 = 0.043 GeV, 

(11) (mjj)p ~ —24.11 GeV + 0.59 m t ; 8 = 0.022 GeV, 

where the subscripts H and P refer to HERWIG and PYTHIA and 8 is the mean squared 
deviation in the fit. Eqs. (10) and (11) correspond to the straight lines in Fig. 4: given the 
slopes of such lines, the difference between HERWIG and PYTHIA A{m Be ) ~ 1.2 GeV 
implies an uncertainty Am t ~ 2 GeV, if one reconstructed nit from a possible measure- 
ment of (mj(). Ref. [23] also pointed out that if one restricted the analysis to the range 
50 GeV < rriBi < 120 GeV and computed truncated moments, the induced uncertainty 
would be reduced to about Am t ~ 1.5 GeV. Since such numbers are comparable with 
the current error on the top mass world average, the conclusion of [23] is therefore that 
it is advisable using the tuned versions of Monte Carlo generators, and ultimately even 
the HERWIG++ code [29], as it fares pretty well with respect to i?-hadron data. 

Before concluding this paper, it is worthwhile discussing the work [30], wherein the top 
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mass is reconstructed as the W + 6-jet combination in the lepton+jets channel, using kx 
and cone (PxCone, infrared-safe) algorithms. It was found out that the spectra obtained 
by using the kr algorithm are mostly affected by initial-state radiation and underlying 
event, whereas the cone algorithm turns out to be sensitive, above all, to final-state 
radiation and hadronization. The recommendation of [30] is that using both algorithms 
will be compelling for the sake of a reliable estimation of the theoretical uncertainty. 

In summary, I discussed a few topics relevant to the top mass reconstruction at hadron 
colliders, taking particular care about the relation with the theoretical mass definitions 
and the uncertainties on the top mass determination. In particular, I presented recent 
work aimed at extracting the pole or MS masses by comparing cross section measurements 
and precise QCD calculations and relating the mass implemented in Monte Carlo event 
generators to the jet mass in the SCET framework, employing i?-evolution. Furthermore, 
progress has been lately achieved in the treatment of bottom-quark fragmentation in 
top decays and Monte Carlo tuning, as well as in the understanding of the sources of 
theoretical error, of both perturbative and non-perturbative origin. 
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